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| According to recent progress in the finite size scaling theory of critical disordered systems, the 

. nature of the phase transition is reflected in the distribution of pseudo-critical temperatures T c (i, L) 

over the ensemble of samples (i) of size L. In this paper, we apply this analysis to the derealization 
transition of an heteropolymeric chain at a selective fluid-fluid interface. The width AT C (L) and 
the shift [Te(oo) — T c a "(L)] are found to decay with the same exponent L _1/V - R , where 1/ur ~ 



0.26. The distribution of pseudo-critical temperatures T c (i,L) is clearly asymmetric, and is well 
fitted by a generalized Gumbel distribution of parameter m ~ 3. We also consider the free energy 
distribution, which can also be fitted by a generalized Gumbel distribution with a temperature 
dependent parameter, of order m ~ 0.7 in the critical region. Finally, the disorder averaged number 
£h ' of contacts with the interface scales at T c like L p with p ~ 0.26 ~ \/vr. 

c ' 

I. INTRODUCTION 

■ 

Hcteropolymers containing both hydrophobic and hydrophilic components are of particular interest in biology 
P|. In a polar solvent, these heteropolymers prefer conformations where the hydrophilic components are in contact 
with the polar solvent, whereas hydrophobic components avoid contacts with the solvent. The behavior of such 
heteropolymers in the presence of an interface separating two selective solvents, one favorable to the hydrophobic 
components and the other to the hydrophilic components, is less obvious, and has been much studied recently 
In the initial work of ref. 0, a model was proposed and studied via Imry-Ma arguments, an analysis of the replica 
Hamiltonian and numerics : it was found that a chain with a symmetric distribution of hydrophobic/hydrophilic 
components is always localized around the interface at any temperature (in the thermodynamic limit), whereas a 
chain with a dissymmetric distribution of hydrophobic/hydrophilic components presents a phase transition separating 
a localized phase at low temperature from a delocalized phase into the most favorable solvent at high temperature. 

^ ■ Experimentally, the presence of copolymers was found to stabilize the interface between the two immiscible solvents ||| , 

since the localization of the heteropolymers at the interface reduces the surface tension. By now, the predictions of ref. 
have been confirmed in the physics community by various approaches includingmolecular dynamics simulations [|| , 
y— I ' Monte Carlo studies 0, ^2 > variational methods for the replica Hamiltonian Q Q , exact bounds for the free-energy 
^0|, and dynamical approach pd| . A real space renormalization group study, based on rare events, was proposed in 
[l3f . Mathematicians have also been interested in this model. The localization at all temperatures for the symmetric 
case was proven in ref. [141 Il5| . In the dissymmetric case, the existence of a transition line in temperature vs 

S dissymmetry plane was proven in ref. 0, 0] . Recently, important progresses were made in various aspects of the 
phase transition [H[ll UliJ. 
. We study here the same model using the Wiseman-Domany approach to finite size scaling in disordered systems 
[22l I23L l24l I2M I2H |27| . The main outcome of this approach can be summarized as follows. To each disordered 
sample (i) of size L, one should first associate a pseudo-critical t emp eratur e TJ i,L), defined for instance in magnetic 
systems as the temperature where the susceptibility is maximum |22Ll24U25ll27| . The disorder averaged pseudo-critical 
J> . temperature T® V (L) = T c (i,L) satisfies 

•i-H , 

X' rrW-^M-L- 1 /^ (i) 

where vr is the correlation length exponent. Eq. Q generalizes the analogous relation for pure systems 

TP ure (L) - T c {oo) - L-V"p (2) 
The nature of the disordered critical point then depends on the width 



AT C (L) = ^/Var[T c (i,L)] = \ T?{i, L) - (T c (i, L) (3) 



of the distribution of pseudo-critical temperatures. When the disorder is irrelevant, the fluctuations of these pseudo- 
critical temperatures obey the scaling of a central limit theorem as in the Harris argument : 

AT C (L) — L~ d / 2 for irrelevant disorder (4) 
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This behavior was first believed to hold in general [22 . 12 1| , but was later shown to be wrong in the case of random 
fixed points : in jUI^, it was argued that at a random critical point, one has instead 

AT C (L) ~ L^ 1 ^ 11 for random critical points (5) 

i.e. the scaling is the same as the L-dependent shift of the averaged pseudo-critical temperature (IJ. The fact that 
these two temperature scales remain the same is then an essential property of random fixed points that leads to the 
lack of self-averaging for observables at criticality [23l |25| . 

In this paper, we apply this finite size scaling analysis to the case of the polymer chain at an interface between two 
selective fluids. The paper is organized as follows. The model is defined in Section|n] The definition of pseudo-critical 
temperatures, and their distribution is presented in Section ITTT1 The distribution of the free energy difference between 
delocalized and localized phases is studied in Section Hvl Finally, in Section we consider the averaged number of 
contacts with the interface. The numerical implementation through the Fixman-Freire procedure is postponed to the 
Appendix. 



II. MODEL AND NOTATIONS 



A. Definition of the model 



The model we consider is defined by the partition function 

Z(L) = ^2 ex P P lasgnzg (6) 

{z Q } V a / 

with inverse temperature (3 — 1/T. The sum in © is over all random walks {z a } of L steps, with increments 
z a+ \ — z a = ±1 and (bound-bound) boundary conditions Z\ = = z^. In eq. JJJJ, the {q a } represent quenched 
random charges, and with our convention, monomers with positive charges prefer to be in the upper fluid (sgnz > 0). 
Furthermore, we have taken <72a+i = for all (a), so that there is no frustration at zero temperature (each monomer 
Z2 a being in its preferred solvent). Even charges q2 a are random and drawn from the Gaussian distribution with mean 
value g (go > 0) 

„/ n 1 l'2a -TO) 2 , . 

From previous work, model 10 is expected to undergo a phase transition at a critical temperature T c ~ 0(A 2 /qo) 
between a localized phase and a delocalized phase in the upper fluid. 



B. Recursion relations for the partition function 



Let Z(a) be the partition function of the partial chain (1, 2, ...a), with monomer (a) sitting at the interface (z a = 0). 
Since the last loop before monomer (a) can be in the upper (+) or in the lower (— ) fluid, we may write 

Z(a) = Z + {a) + Z-(a) (8) 

The partition functions Z (a) satisfy the recursions 

Z ± (a)= ^ exp[±l3 1i] N{a!,a)[Z + {a') + Z-{a')] (9) 

a'=l,3,.. V i=a' + l J 

where J\f(a',a) represents the number of random walks of (a — a') steps going from z a > — to z a — 0, without 
touching the interface at z = in between. For the critical properties of the derealization transition, the only 
important property of J\f(a',a) is its behavior at large separation 

Af(a',a) ^ a—?—— (10) 
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where a is a constant of order 1. In our numerical study, we have chosen to take the following form 

2a — a' 



J\f(a' , a) = a 



(a - a' - l) 3 /2 



(11) 



for all (a — a'). If one forgets the boundary conditions at (1, a), the partition function characterizing the delocalized 
phase in the (+) solvent (qo > 0) would simply be for each sequence of charges 



Z deloc (a) = 2 a ~V y ( Q ) with V(a) = J2 ® 

a' =2 

It is convenient to introduce the ratio 

where Z^(a) = Z ± (a)/Z deloc (a) satisfy the recursion relations (Eq. [jJJ 



(12) 



(13) 



^ (a -a' -l) 3 / 2 

a'=l,3... v ' 



Z+(a') + Z-(a') 



,3 

OL-1 



a' = l,3,. 



(a - a' - 1)3/2 



Z+{a') + Z-(a') 



(14) 



The numerical iteration of Eqs l|14|) takes a CPU time of order Q(L 2 ) for a chain of length L, but it is possible 
to obtain a CPU time of order O(L) via a Fixman-Freire scheme [28| . Originally proposed in the context of the 
Poland-Scheraga model of DNA denaturation, this scheme can be adapted to the interface model, as explained in the 
Appendix. 



C. Numerical parameters 

We have taken qo = 1 and A = 10 in the Gaussian distibution of eq Q. The data we present correspond to various 
sizes L with corresponding numbers n s (L) of independent samples. Unless otherwise stated, we have considered 
the following sizes going from L = 2 • 10 3 to L = 2048 • 10 3 , with respectively n s {L = 16 • 10 3 ) = 480 • 10 3 to 
n s (L = 2048 • 10 3 ) = 15 • 10 3 . More precisely, we consider 

16, 32, 64, 128, 256, 512, 1024, 2048 (15) 
480,240,120,60,30,15,30,15 (16) 



L 
1000 

n s {L) 

1000 



III. DISTRIBUTION OF PSEUDO CRITICAL TEMPERATURES 



In a previous paper 29] , we have discussed two possible definitions for the pseudo-critical temperatures for Poland- 
Scheraga models. Here we use the definition based on the free energy, as we now explain. 



A. Definition of a sample dependent T c (i,L) 

For each sample (i) of length L, we may define a pseudo-critical temperature T c (* ; L) as the temperature where the 
free energy F« (L, T) = -TbxZ^(L,T) reaches the delocalized value F^ loc (L, T) = -T(i-l)ln2-V(i) (Eq. EJl, 
i.e. T c {i, L) is the solution of the equation 



(L, T) + T(L - 1) In 2 + V(L) = 



(17) 
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As explained in details in our previous study [23, this definition of pseudo-critical critical temperatures, together 
with bound-bound boundary conditions, introduces logarithmic correction in the convergence towards T c (oo). Eq. 
JU is accordingly replaced by 



2r(L)-T c (oo)~ 



(18) 



B. Scaling form of the probability distribution 

Our data for the distribution of pseudo-critical temperatures (see Figure ^) follow the scaling form 



Pl(T c (i,L)) 



AT C (L) 



g\x = 



T c (i ,L)-T™(L) 
AT C {L) 



(19) 



where the scaling distribution g(x) (normalized with < x >= and < x 2 >= 1) is well fitted by the one parameter 
generalized Gumbel distribution G m (— x) 



g(x) = G m (-x) 



r(m) 



(.--)• 



(20) 



with parameter m ~ 3. On Fig. ^b), we show the fit of our data with G^{— x), and plot for comparison the usual 
Gumbel (m = 1) and Gaussian (m — > 00) distributions. 





FIG. 1: (a) Distribution of critical temperatures for (L/10 3 = 32, 64, 128, 256, 512, 1024, 2048) (b) Same data after the rescaling 
of eq. (1191 : fit with the generalized Gumbel distribution G?3(— x) (thick line). For comparison, we also show the usual Gumbel 
m = 1 (thin line) and Gaussian m = 00 (dashed line) distributions. 



C. Scaling properties of the shift and of the width 

We now consider the behavior of the width AT C (L) and of the average T^ V (L) of the distribution l|19|) . as L varies. 
For the width AT C (L), we have fitted our data with the following power law (see Fig|2Ja)) 

. x 0.26 



AT C (L) ~ 86 (^-J (21) 
For the average T c a "(L) (see Figl^b)), our data are well fitted by the generalized form of eq. ljT%|l 

IT (L) * 83.8 - 246 ( MV 289 ) ) 26 (22) 
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Log L L 

FIG. 2: (a) Log-Log plot of the variance Var[T c (i,L)] = (AT C (L)) 2 for L/10 3 = 32,64,128,256,512,1024,2048 (O) with the 
fit of equation 12 II (b) The average T^ V (L) for the same sizes, with the fit of equation 1221 . 



The value T c (oo) ~ 83.8 = 0.838-^- is in agreement with the numerical estimate of ref .20j. The result for the exponent 
in eqs l|21|) and l|22|) indicate that the transition can be described as a random critical point with a single correlation 
exponent 

— ~ 0.26 (23) 

V R 

A similar value has been observed in numerical simulations by the authors of ref. |20| (private commnunication 
Our data rule out the possibility of an infinite order transition based on rare negatively charged sequences 
The present results suggest that the excursions in the unfavorable fluid that are important for the transition, are of 
finite length. 

IV. FREE ENERGY DISTRIBUTION 

Since the distribution of pseudo critical temperature is clearly linked to the statistical properties of the free energy, 
we have studied the probability distribution of the free energy difference F^(L,T) defined by (see eq. 117(1 

F®(L,T) = F®(L,T) - F^ loc (L,T) = F®(L, T) + T(L - 1) In 2 + V(L) (24) 

Defining the average F^ av \L,T) and the width AF(L,T), we write the scaled probability distribution as 

P(FW( L ,D) * - 1 h T L = ^(L,T)-F^{L,T)\ 
AF(L,T) \ AF(L,T) J 

where the function Iit(x) is independent of L, but depends on the temperature T. We show in Figure the function 
}it{x) for T = 40 (deep in the localized phase), T = 82 (critical region) and T = 130 (deep in the delocalized phase). 
It turns out that a good one parameter fit is again obtained with the generalized Gumbel distribution 

h T {x) = G m {x)^^(e*- e:c ) m (26) 

where the effective parameter m depends on the temperature T. We now characterize the behaviour at low and high 
temperatures, as well as in the critical region. 

A. The localized phase 

At T <C T c , the sample dependent free energy can be fitted by 

F i (L,T)=Lf 0O (T)+a T VLu i + 0(l) (27) 
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FIG. 3: (a) Fit of the rescaled free energy distribution kr(x) Cea. ^ 1 at T = 82 for L/10 3 = 32,64,128,256,512 by the 
generalized Gumbel distribution G m {x) with m — 0.66 (thick line). The usual Gumbel distribution (m = 1) is also shown for 
comparison (dashed line) (b) The function h T (x) for L/10 3 = 32 at T = 40 (thick solid line), T = 82 ~ T c (dashed line) and 
T = 130 (thin solid line). Fitting with the generalized Gumbel distribution G m {x) leads to a non monotonous variation of m 
w.r.t the temperature T. The minimal value of m is obtained in the critical region. 

where Ui is a Gaussian random variable (which corresponds to the Gumbel parameter m — > oo). Eq. (|27[1 is indeed 
expected from the existence of an extensive number of contacts with the interface: fluctuations are then of order \fL 
and Gaussian distributed. 

B. The delocalized phase 

At T ^> T c the sample dependent free energy Fi(L,T) can be fitted by 

Fi(L,T) = ~T La- +vi (28) 

l Lq 

where Iq is a constant, and Wj is a random variable distributed with G m (x) and to ~ 1.3 (sec cq. (I26[) . The first 
term on the r.h.s of ea.(|28fl comes from the bound-bound boundary conditions, while the random term of order 0(1) 
reflects rare contacts with the interface. 

C. In the critical region 

In the critical region, the average free energy F^ av '(L, T) behaves as 

FW(L,T) = LfoeiT) + /i(T) Ln^ (29) 

where fi is a coefficient which weakly depends on temperature in the critical region. The thermodynamic free energy 
density foo(T) is then found to go to zero between T = 82 and T = 83, in good agreement with the value T c ~ 83.8 
obtained above. Concerning the rescaled distribution hxix), we observe that the parameter m remains at a minimum 
value of order 0.7 in a large temperature range 75 < T < 85. This large crossover region is clearly linked to the small 
value of 1/vr ~ 0.26. 

V. CONTACTS 

We have also considered the disorder averaged number of contacts N with the interface in the critical region, and 
we numerically obtain 

L 

N (T C , L) = J2 <S(z a ,0) > Ta ~ L p , p 0.26 (30) 

a=l 
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It is interesting to compare this result with the behavior of the internal energy difference Eq(T,L) between the 
localized and delocalized phases. A simple scaling argument based on the hyperscaling relation a — 2 — vr yields 
the following behavior for the energy Eq(T Ci L) ~ L X I VR at criticality. The fact that the exponent p is close to the 
exponent l/V/j points towards the same scaling for the energy and number of excursions in the unfavorable fluid : 
this means that, at criticality, the excursions in the lower unfavorable fluid are of finite length. 

Finally, our data on the number of contacts Nq(T, L) yield an estimate of the critical temperature 82.5 < T c < 83, 
which is again compatible with the above mentionned values. 



VI. CONCLUSION 



We have studied the probability distribution finctions (PDF) of the pseudo-critical temperatures T c (i, L) and of the 
free energy Fi(L,T) for the localization transition of an heteropolymer at a selective fluid-fluid interface. We have 
obtained an estimate of the critical temperature T c ~ 0.83^-, in agreement with J2jj. Both the width AT C (L) and 

the shift [T c (oo) — T£ V (L)] are found to decay with the same exponent L~ X / UR where 1/vr ~ 0.26. These results show 
that the transition is not due to rare events, such as large excursions in the unfavorable fluid |l3j |. 

Concerning the form of the PDF's, we find here asymmetric distributions, which can fitted with generalized Gumbel 
distributions G m (x). This is in contrast with the wetting and Poland-Scheraga transitions [29| . where these PDF's were 
found to be Gaussian. Generalized Gumbel distributions have been recently used in various contexts 
As discussed in ref. |3l| . it is empirically known that PDF's with the same first four moments approximately coincide 
over the range of a few standard deviation which is precisely the range of numerical -or experimental- data. So 
generalized Gumbel PDF's, with arbitrary m, should probably not be considered more than a a convenient one 
parameter fit (for instance, Tracy- Widom distributions defined in terms of Airy functions 36], can be also fitted very 
well by generalized Gumbel distributions in the range of interest). 
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APPENDIX A: FIXMAN-FREIRE SCHEME FOR THE INTERFACE MODEL 



The Fixman-Freire scheme |28|. invented in the context of DNA denaturation, allows one to compute melting 
curves of DNA sequences in a CPU time of order 0(L) instead of 0(L 2 ). Here we briefly explain how to adapt the 
Fixman-Freire method to the interface model. 

The first step consists in replacing the factor l/(a — a' — l) 3 / 2 appearing in the recursion relations (Eq. I14|) by 

l^J^e-W— (Al) 

(a — a' — 1 ) A I Z r—' 

where the number / of couples {(ii,bi) depends on the desired accuracy. The parameters (di,bi) are determined by a 
set of non-linear equations. This procedure has been tested on DNA chains of length up to L = 10 6 base pairs [3?ll3q |. 
and the choice / = 15 gives an accuracy better than 0.3% and we have adopted this value. The table of parameters 
(a,i,bi) we have used can be found in our previous work on the wetting transition |30j |. 
Let us now introduce the following notations 

RHa) (A2) 

-L ] =R+( a )+R-{a) (A3) 

S{a) Z(a) 
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as well as the coefficients (if (a) defined in terms of the Fixman-Freire coefficients (bj) (Eq. I Al|) where j = 1, 2, = 15 



a— 4 

/ut(a-2) = e -M Q - 4 ) ^ e b ^'S(a-2)S(a-4)...S(c/ + 2) 

a' = l,3,.. 

a— 4 



M7(«-2) = 



_ „-6 J (a-4)-2/3y(a-4) 



^ e ^a'+2/3y( Q ')5( a _ 2 )S( a - 4)...S(a' + 2) 



a'=l,3,.. 

It is then easy to obtain that for any j, these coefficients satisfy the simple recursions 



4(a) =S(a)[l + e-^4(a-2)] 



IMj(a) = S(a) [1 + e-^- 2 ^'-'iMj(a - 2)] 
In terms of these coefficients, the recursions (Eq. I14|l become for the ratios i? ± (a) (Eq. IA3(1 



R+( 



a) = a 



l + ^ 0i e- 3 ^/i+(a-2) 



The recursions (|A7IA8|I can now be iterated in a CPU time of order O(L) 



1 + e -W**- 3 J2a j e- 3b * f ij(a-2) 

3 = 1 



(A4) 
(A5) 



(A6) 
(A7) 



(A8) 
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